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We study the formation of shocks on the surface of a granular material draining through an orifice 
at the bottom of a quasi-two dimensional silo. At high flow rates, the surface is observed to deviate 
strongly from a smooth linear inclined profile giving way to a sharp discontinuity in the height of the 
surface near the bottom of the incline, the typical response of a choking flow such as encountered in 
a hydraulic jump in a Newtonian fluid like water. We present experimental results that characterize 
the conditions for the existence of such a jump, describe its structure and give an explanation for 
its occurrence. 



The flow of granular materials presents a multitude 
of experimental and theoretical challenges in many-body 
physics, as manifested in the unusual collective behavior 
of a large assembly of macroscopic grains that interact 
with each other through collisions and friction. The phe- 
nomena of interest in quasi-static motion range from jam- 
ming transitions, the formation and evolution of shear 
bands, and the transmission of forces and acoustic signals 
in granular packing . In the dynamical regime asso- 
ciated with technological applications such as silo flows, 
other questions of interest include the origin and mainte- 
nance of fluctuations and their relation to mean flows in 
the bulk, and the nucleation and evolution of avalanche 
like behavior in the vicinity of the surface ■ 

A system of long-standing technological interest in this 
last situation is associated with flow in a silo. Build- 
ing on an understanding of simple properties such as the 
flow rate as a function of the draining orifice size and 
the shape of the silo g , recent work has focused on the 
details such as the form of the velocity profile of the bulk 
flowing regions. In we studied the flow in the vicinity 
of the orifice of a silo. The predominant flow is restricted 
to a parabolic region centered at the origin and may be in- 
terpreted in terms of simple diffusion models (||| ■ Here, 
we complement those studies by considering the dynam- 
ical evolution of the form of the free-surface in a quasi 
two-dimensional silo flow. 

Our system consists of a flat-bottomed silo that is 
quasi-two-dimensional: 89 cm wide, 45 cm high and 2.54 
cm deep, the last being between 25-1000 times the typi- 
cal grain size. The orifice itself has a square cross-section 
of width equal to that of the silo, with an attached valve 
to control the flow rate Q (see Fig. |l|). We use spherical 
glass grains with diameters 1 mm, 500 /im, 100 /im, and 
30 /xm to study the dependence of the phenomenon on 
grain size. 

The surface and the bulk of flow is visualized through 
the glass side walls of the silo using two types of CCD 
cameras. The Kodak ES 1.0 digital camera with a res- 
olution of 1000 pixels x 1000 pixels with a maximum 
frame rate of 30 per second was used to measure the sur- 



face profile and to determine the location of the static 
and mobile regions inside the silo. A second high speed 
Kodak camera with a maximum frame rate of 1000 per 
second and resolution of 256 x 240 pixels was used to mea- 
sure the velocity of the particles. By focusing in small 
regions where a particle corresponds to 3-4 pixels, we are 
able to identify and track individual particles over con- 
secutive frames with an accuracy of < 0.1 pixels and thus 
determine the velocity. 

When the orifice is opened, bulk flow is rapidly estab- 
lished and the free surface initially subsides as an inverted 
Gaussian and then quickly develops in to a V shape. The 
angle of inclination of the surface approximately corre- 
sponds to the angle of repose of the glass beads. For low 
flow rates, the free-surface has a constant slope, as shown 
in Fig. |2|(a). However as the flow rate is increased past a 
threshold by opening the valve, the inclined free-surface 
is no longer linear but shows a region of abrupt height 
transition right above the orifice, as shown in Fig. ^|(b). 

This unusual surface feature is similar to a hydraulic 
jump in a Newtonian fluid observed when a stream of liq- 
uid encounters an obstacle (such as observed when water 
impinges onto a flat plate) . Past work JlO|— has shown 
that shock-like features similar to hydraulic jumps are 
observable in granular flows in the presence of an exter- 
nally introduced physical obstacle such as a vertical plate 
in a channel flow. In these studies it was noted that the 
the ratio of inertial force to the gravitational force in the 
system given by the Froude number exceeds one, in anal- 
ogy to the hydraulic jump observed in a fluid. Although, 
in granular flow it is less obvious that this is important 
as in fluids where it corresponds to the flow velocity ex- 
ceeding the wave speed for shallow depths. This fact was 
also noted earlier [nj. 

In the case of the silo flow treated here, the jump is 
spontaneously generated by the flowing grains which cre- 
ates its own "obstacle" by choking right above the orifice. 
Moreover, unlike in a Newtonian fluid or in the channel 
flows of granular materials where the entire fluid is in 
motion, here the grains exist in different "phases" in dif- 
ferent regions; away from the orifice and the free-surface, 
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they barely move at all, in the neighborhood of the ori- 
fice they move slowly, near the free-surface they flow like 
liquids, and veritably boil into a gas in the vicinity of the 
free-surface shock [Fig. ||(c)]. To understand the forma- 
tion of the free-surface shock it is important to under- 
stand the nature of the granular flow below the surface. 
We locate the surface and the boundary between the mo- 
bile and static regions as the silo empties using the high 
resolution digital camera. Two temporally separated im- 
ages are subtracted to determine the boundary between 
static and mobile regions. This boundary along with the 
surface is shown for one-half of the silo in Fig. || as the 
granular matter discharges. 

Since all the grains must eventually exit through the 
orifice, the vertical velocity v{x,y) is largest there, and 
the horizontal velocity u(x, y) is smallest. However, as 
one moves away in the vertical direction v gets smaller, 
while u correspondingly becomes larger, and the extent 
of the flowing region becomes larger. As voids diffuse up- 
wards from the orifice, they explore a larger and larger 
distance in the horizontal direction, and the bulk flow 
seems to be driven by diffusive motion. In fact, the flow 
in the silo close to the orifice can be described remarkably 
well in terms of an invariant boundary between flowing 
and non-flowing regions of the form y ~ x 2 j d with d be- 
ing the grain diameter Q , consistent with the kinematic 
model for silo flows ||. 

The flow near the free surface outside this parabolic re- 
gion is confined to a shallow layer where particles move 
essentially parallel to the surface. However above the 
orifice, the flow near the free surface changes from pre- 
dominantly along the surface to a bulk flow towards the 
orifice. The size of this intermediate region is a function 
of the height h of the free surface above the orifice; this 
region is very large at early times, and gets progressively 
smaller. We parameterize this intermediate region us- 
ing the distance I from the bottom of the surface to the 
point where the flow can be described by the parabola 
(dashed line in Fig. ||). The inset to Fig. ^ shows I as 
a function of h. Using the high speed camera, we mea- 
sured the velocity of particles v s (s, z) in the neighbor- 
hood of four equidistant locations along the free surface, 
with Si,i — 1 — 4 indicated in Fig. |^ as a function of the 
depth z perpendicular to the surface. In Fig. |], we plot 
v s (si,z) vs. z on a log-linear scale. The velocity decays 
rapidly for the shallow flowing regions (si,sa) and is de- 
scribed by v s (s, z) = vq exp(— a\z/d — a2(z/d) 2 ) where 
ct\ and a2 are dimensionless quantities that depend on 
the location of the surface. In the intermediate region 
(S3, S4), the velocity profile can be described by a single 
exponential v s (s, z) = vq exp(— ctiz/d). Fig. ^ also shows 
that the particles quickly reach a steady velocity charac- 
terized by the balance between gravitational forcing and 
inelastic collisions p3[ . 

We followed the free-surface as a function of time using 
the high resolution camera (Fig. m while simultaneously 



measuring vq (Fig. ||) for the 100 /im glass beads. In 
Fig. ||, we see that vq increases as the silo empties. For 
low Q, the surface always remains V shaped throughout 
the draining process. For larger Q the V shaped surface 
changes into a U shaped surface eventually leading to 
the formation of a granular jump. For 100 /mi beads the 
shock appears at a height of about 10 cm for a Q of 250 
g s _1 , and at a height of about 6 cm for a Q of 170 g s _1 . 
A shock is not observed if the Q is decreased further by 
using the valve. Thus the shock appears only when the 
surface velocity exceeds a critical value and occurs at a 
higher height for higher Q (dashed line in Fig. |6|). 

We also used beads of various sizes to study its effect on 
the appearance of the shock. It has been proposed that 
the flow rate Q oc (W — ad) 5 / 2 , where W is the width of 
the silo, and a is a constant between and 1 depending 
on the geometry of the grain with size d Q. However, we 
observe that the size of particles has a greater influence 
on the flow rate than accounted for by this formula. The 
maximum Q for 30 //m, 100 //m, 500 /mi and 1 mm beads 
was 260 g s _1 , 250 g s _1 , 230 g s" 1 and 200 g s -1 respec- 
tively for the same (2.54 cm x 2.54 cm) orifice. In case 
of 30 /im beads, the shock appears when the surface is 
approximately 30 cm above the orifice as opposed to the 
10 cm for 100 /im. In case of 1 mm beads, a clear shock 
does not develop although a stagnation zone at the base 
of the incline surface is observed. We plot the critical 
height h c versus the flow rate for various d in Fig. [?| to 
characterize the transition of the surface from V-shape 
to that where a jump is observed. We note that the ratio 
of the drag force to weight of the particles mg, is less 
than 0.1% for 1 mm particles and less than 1% for 100 
/im and therefore negligible. Therefore we believe that it 
is not the cause of the increase in height. The drag force 
for 30 /im particles is about 10% and may be important. 
The previous experiments ffjj for 2.54 cm and 1.27 cm 
wide silos, showed that the wall friction does not signifi- 
cantly add to the drag force experienced by the particles 
as measured by the angle of repose and the shape of the 
parabolic flow region and therefore is negligible. 

To rationalize these observations, we start by consider- 
ing the increase in vq which occurs for essentially geomet- 
ric reasons. The flow rate Q is proportional to the grain 
flux past the dashed lines shown in Fig. |j| and scales 
a,s W x J v s (s,z)dz, where W is the width of the silo. 
Using the observed exponential dependence of v s (s,z), 
Q ~ vq x ho. As the silo empties, the bulk flow is con- 
fined to a smaller and narrower region (see Fig. ||). We 
have observed Q by measuring the mass of the granular 
matter flowing out of the orifice in a fixed interval of time 
and find that Q is constant to within 5% O] . Since Q is 
constant, vq must increase with time in the intermediate 
region whose extent simultaneously decreases. Therefore 
the transition from surface to bulk flow becomes more 
and more abrupt. This increases the frequency of inelas- 
tic collisions in the neighborhood of the valley leading to 
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a region of lower mean velocities (see Fig. p](c)). The de- 
crease in the horizontal component of the velocity as the 
grains approach the center of the silo leads to the condi- 
tions for a kinematic shock as faster grains pile up behind 
(and atop) the slower grains. Owing to the symmetry of 
the flow from either side, the free surface responds by first 
flattening along its axis of symmetry (see Fig. ||). As vq 
exceeds a critical value, an abrupt increase in the height 
of the free surface is the only recourse for the rapidly 
flowing surface grains which form a granular jump. 

In case of collisional granular flows, theoretical ap- 
proaches have been developed first by M| and jl6) using 
mass-momemtum-energy-conservation to derive contin- 
uum equations. These approachs and subsequent im- 
provements have been also reviewed by [fL7f . However, 
applying these continuum approachs, especially in case 
of free boundary conditions as we have on one side of the 
flow and enduring contacts at the bottom of the surface 
flow, is extremely difficult. Therefore in the following we 
attempt scaling arguments to explain the phenomena. 

To quantify the criterion for the formation of the jump, 
we consider the analogous case of a hydraulic jump in 
a fluid which occurs when the velocity of flow exceeds 
the propagation speed of a surface perturbation. Di- 
mensional arguments lead to an estimate of the speed 
of propagation of surface perturbations as v p ~ \J gD, 
where D corresponds to the depth of the flow [lq ]. In 
our case D ~ 20 x d. Therefore v p ~ 14 cm s" 1 for 100 
fj,ra beads, in qualitative agreement with the surface ve- 
locity at which the shock is observed to occur. Finally, 
because v p increases with d, this is also consistent with 
the occurrence of the shock at a much higher height for 
30 /urn particles and the lack of a shock in case of the 
larger 1 mm beads. 

In conclusion, we have shown that a dynamic jamming 
event in silo flows can lead to the formation of free sur- 
face shocks that are similar to hydraulic jumps in simple 
fluids whose flow is obstructed. An interesting feature 
of the granular jump is that the obstruction is generated 
spontaneously by the flowing grains. We have given a 
kinematic explanation for our observations, but in the 
absence of an elaborate theory that describes granular 
flows, a quantitative explanation remains an open ques- 
tion. 
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FIG. 3. The surface and the boundary between flowing 
and static regions at three different times. The flow region 
is parabolic, and stay parabolic near the orifice during dis- 
charge. The orifice is centered at the origin (0, 0). The dashed 
lines show the size of the intermediate region. The points sf , 
s2, s3, and s4 mark the locations where the velocity is mea- 
sured and plotted in Fig. d = 1 mm. Inset: The length I 
that measures the size of the intermediate region as a func- 
tion of height h of the surface above the orifice is observed to 
decrease as the silo empties. 



2.54 cm 

FIG. 1. The experimental apparatus is a quasi- two dimen- 
sional silo with glass side walls. The dashed curve indicates an 
example of the boundary between rapidly flowing and static 
regions inside the silo. The dotted square box close to the ori- 
fice shows the approximate area corresponding to the images 
shown in Fig. |^(a) and Fig. ^(b) . 



FIG. 2. (a) Image of granular matter draining from an ori- 
fice inside a quasi- two dimensional silo as a function of time. 
The surface is linear in this case. Q ~ 50 g s"\ d = 100/im. 
The images correspond to a 7 cm x 7 cm region close to the 
orifice, (b) The surface is observed to show an increase in 
height similar to a hydraulic jump when a rapidly flowing 
fluid in a channel encounters an obstacle. Q = 200 g s _1 . (c) 
Closeup image ( x 3) of the shock as it is forming. The dura- 
tion of the exposure is 4 ms and therefore the fast moving par- 
ticles appear as streaks. By comparison, the particles inside 
the jump appear almost frozen. A movie of the jump can b e 
viewed at jittp : / /physics . clarku . edu/ ~akudrolli / shock. html , 
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FIG. 4. The velocity v a (s,z) as a function of depth z at 
four positions si — S4 along the surface (indicated in Fig. |^). 
d = 1 mm. The s-axis is chosen parallel to the surface and 
the z axis is chosen perpendicular to the surface and pointing 
inwards. The velocities are observed to decay approximately 
exponentially from the surface. The length over which the 
velocity decays ho is approximately 10-20 grain diameters. 
Note also that the particles have reached an approximately 
constant velocity v s (s, 0) = Vq along the surface. 




FIG. 5. The surface profile of the granular matter draining 
from a silo as a function of time showing the development of 
the jump. (Q ~ 170 gs , d— lOOpm). The origin is located 
at the center of the orifice. The first four surface profiles are 
separated by 5 second time intervals, to show the features of 
the jump, the last five are separated by 3 second time inter- 
vals. The dashed lines are guides to the eye to compare the 
surface profile to that at low Q. The surface flattens and then 
a shock is observed to develop. 
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FIG. 6. The velocity of grains vo at the surface as a func- 
tion of height above the orifice for d = 100/im and d — 1 mm 
glass beads. The velocities are measured at a position just 
before the location where the jump occurs (S4 in Fig. |^). The 
dashed line corresponds to the estimate of the critical velocity 
of 14 cm s _1 for 100 micron particles. The critical velocity 
for 1 mm particles is approximately 44 cm s~ . 
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FIG. 7. A plot of the critical height h c where the jump first 
appears as a function of the flow rate Q for glass particles of 
various diameters d. A clear jump is not observed for 1 mm 
particles athough a small stagnant region is observed at the 
fastest flow rate. 
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